Fermi-surface reconstruction in a smectic phase of a high temperature superconductor 
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It is shown that, in the presence of a moderately strong d symmetry breaking (which could be 
produced either by lattice orthorhombicity or the presence of an electron nematic phase), a weak, 
period 4, unidirectional charge density wave ("charge stripe") order can reconstruct the Fermi surface 
of a typical hole-doped cuprate to produce a small electron pocket. This form of charge density 
wave order is consistent with that adduced from recent high field NMR experiments in YBCO. The 
Fermi pocket has an area and effective mass which is a rough caricature of those seen in recent high 
field quantum oscillation experiments. 

PACS numbers: 



There is clear evidence that, in addition to supercon- 
ductivity, other broken symmetry electronic phases occur 
in some circumstances in the cuprate high T c supercon- 
ductors. Discovering and characterizing these orders is 
an integral part of cuprate research. 

A high level of excitement has accompanied the 
discovery [lj of quantum oscillations in high quality "un- 
derdoped" crystals of YBCO when superconductivity is 
quenched by high magnetic fields 0-0] . What is appeal- 
ing about these experiments is that they potentially re- 
veal the properties of ordered groundstate phases which 
compete with superconductivity. Specifically, these ex- 
periments unambiguously establish the existence of small 
Fermi pockets with well defined gapless Fermionic quasi- 
particle excitations. Significantly this "fermiology" is 
quite different from that of highly overdoped materials 
as also inferred from quantum oscillations ^ . Although 
there is some uncertainty concerning the number of Fermi 
pockets in underdoped YBCO, there is growing consensus 
that there is at least one "electron-pocket" which encloses 
roughly 2% of the Brillouin Zone (BZ). While there are 
some reasons to worry Q that this pocket may reflect a 
material specific band which is irrelevant to the physics 
of the copper-oxide plane, it is widely considered more 
likely that the electron pocket is the result of a Fermi 
surface reconstruction brought about by some form of 
translation symmetry breaking density wave order. 

Theoretical calculations have shown that basic features 
of the above fermiology can be accounted for under the 
assumption that the system manifests d-density wave 
(dDW) order 0, spin-spiral order[Ioj], or spin-stripefll| 
order [12J. In contrast, it has been argued [12| that unidi- 
rectional CDW (or charge stripe) order cannot produce 
the requisite electron pocket. 

In this context, it is si gnif icant that a recent high field 
copper NMR experiment [13[ performed on a YBCO sam- 
ple which does exhibit quantum oscillations has lead to 
the following conclusions: 

1) There is no effective magnetic field on the Cu sites - 



this rules out any form of SDW order, although, for spe- 
cial symmetry reasons, it does not rule out dDW order. 
2) There is a charge density modulation most probably 
associated with a site centered unidirectional CDW with 
period equal to 4 lattice constants Moreover, the 

amplitude of charge modulation is estimated to be about 
0.03 ± O.Ole per planar unit cell. 

Very similar forms of charge order have been observed 
directly [ll, 15- IH with neutron and X-ray scattering for 
some time in the 214 family of high temperature super- 
conductors, especially ones with generally reduced su- 
perconducting transition temperatures. Such order is of- 
ten referred to as "charge stripes." Other studies have 
revealed a strong enhancement of spin-stripe order in 
LSCO when superconductivity is suppressed by a moder- 
ate magnetic field[3]. In addition to the above, evidence 
of a tendency to formation of charge density wave order 
with period near 4 lattice constants has been adduced 
from STM studies of other underdoped cuprates[lll [20l — 
[23[, and of field induced CDW correlations in vortex 
cores [2i| . Given these results, it is important to revisit 
the issue of whether a pure uni-directional CDW order 
can produce the requisite electron pocket in YBCO. 

In the following we demonstrate that charge stripes can 
produce an electron pocket with the indicated character- 
istics so long as, in addition to the CDW order, the elec- 
tronic structure also breaks C4 rotational symmetry suf- 
ficiently strongly. Of course, unidirectional order, itself, 
breaks this symmetry, as does the orthorhombic structure 
of YBCO. However, if these explicit symmetry breaking 
effects are too weak, as existing calculations suggest, then 
the strong C4 symmetry breaking could also reflect the 
existence of an intrinsic electronic tendency to form an 
"electron nematic" state[25j, for which there is consider- 
able independent experimental evidence II, 23, 2?j|-3tl|. 

In order to convey our message clearly let us consider 
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the following simple tight-binding Hamiltonian. 



t> [ct 
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- £ [2y cos(7rX/2) + 2F 2 cos(TrY) + M ] cj^c^. (1) 

_R,<7 

For ^jv = Vb = V2 = the model is a caricature of 
the low energy bandstructure associated with the copper 
oxide planes inferred from ARPES studies - fits to the 
Fermi surface shape in Bi2201 [3lJ yield values of t'/t 
between —0.2 and —0.3, and in LSCO [32| between —0.3 
and —0.4. Often these fits also invoke a third-neighbor 
hopping matrix, t"/t between 0.1 and 0.2. For the ex- 
plicit calculations we have carried out, we have taken 
t'/t = -0.3 or -0.4 while t" /t = 0, although within 
some broad range, the results are not highly sensitive to 
the choice of band parameters. 

Because YBCO cleaves between highly dissimilar 
planes, its surface can potentially have a large effect 
on the surface electronic structure accessible to ARPES 
studies. The analysis of the ARPES spectrum is further 
complicated by the (possibly substantial) bilayer split- 
ting of the bands expected in YBCO, and the presence 
in some form of chain-related bands. Despite all this, 
most ARPES studies [33Tl35| have inferred a rather sim- 
ilar "normal-state" Fermi surface structure for YBCO 
as those seen in other cuprates with an inferred value 
of t'/t w -0.4. However, a recent ARPES study[H] of 
YBCO films has observed a Fermi surface shape which 
is rather different from other ARPES results and, corre- 
spondingly, from our assumed model. 

The remaining parameters in H represent various bro- 
ken symmetry. </>aj- is a measure of the degree of asymme- 
try upon 90° rotation in the absence of the density wave 
order, which reflects both the roughly temperature inde- 




FIG. 1: The schematic representation of the Fermi surface 
with (solid line) and without (dashed line) nematic ordering. 
In (a), t'/t = -0.3 and <f> N = 0.074; in (b), t'/t = -0.4 and 
4>n = 0.143 in (a) and (b) respectively. In both cases, the 
doping is x — 1/8. 



pendent effects of the crystalline anisotropy and a poten- 
tially larger and more strongly temperature dependent 
contribution from any intrinsic tendency toward electron 
nematic order [rjj . The term proportional to Vq is the 
fundamental Fourier component of the effective potential 
representing translation symmetry breaking correspond- 
ing to a period-4 uni-directional CDW; the choice of the 
phase of the cosine corresponds to site-centered charge- 
stripe order (as indicated by the NMR results). The term 
proportional V% is generally present as the only allowed 
harmonic of the fundamental period in a site-centered 
charge-striped state; in the special case where we are 
dealing with the Ortho II phase of YBCO, there is a 
further contribution to V2 (which is unrelated to density- 
wave formation) which comes from the existence of alter- 
nating oxygen-full and oxygen-empty chains in the Ortho 
II crystal structure. The chemical potential \x is adjusted 
so that the density of doped holes per planar Cu is in the 
neighborhood of 10%. 

In Fig. [1] we show the Fermi surfaces associated 
with an isotropic normal phase, 4>n — Vo = V2 = 
(dashed line), and with an anisotropic uniform phase, 
0jv + 0,V = V 2 = (solid line) where <f> N = 0.074 
and 0.143 in Figs, la and lb, respectively. Note that we 
have taken <f>^ larger than the critical value <j) c (which are 
4> c = 0.072 and 0.142 for the choice of band parameters 
in Figs. 1(a) and 1(b), respectively) at which point the 
Fermi surface topology changes from closed (and hole- 
like) to open. This sort of Fermi surface reconstruction 
has been invoked previously [38j to account (in the con- 
text of a weak-coupling RPA treatment) for the nematic 
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FIG. 2: Schematic representation (thick black lines) of the 
folded Fermi surface with nematic ordering and period-4 
charge stripe order and x = 1/8. In (a), t'/t = —0.3, 
4>n = 0.074, V = 0.15i, and V> = O.lOt, while in (b), 
t'/t = -0.4, (j> N = 0.143, Vq = O.llt and V 2 = 0.09t, and the 
electron pocket encloses roughly 1% of the unreconstructed 
BZ. The parameters have been chosen for graphical clarity. 
When comparing with experiment we have taken Vo — 0.02t, 
and V2 = O.Olt, so that the electron pockets are roughly 2% 
of the BZ, but then the gaps between Fermi surface segments 
are quite small. The thin (red) lines are the zone-folded Fermi 
surfaces in the limit of vanishing stripe order, Vo = V2 = 0, 
and the (blue) shaded line is the unfolded Fermi surface. 
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order detected [27| in inelastic neutron scattering studies 
of somewhat more underdoped YBCO than the mate- 
rial which exhibits quantum oscillations. Two important 
wave- vectors which characterized the open Fermi surface 
are q and q' ', the minimum and maximum spanning vec- 
tors of the Fermi surfaces in Fig. 1. Clearly, q and q' 
depend on the magnitude of <^jy, the doping level, and 
other band-structure parameters. 

In Figs. 2(a) and 2(b) we show the reconstruction of 
the Fermi surface produced by a non-zero CDW poten- 
tial, Vo,2 > 0. We report the parameters used to con- 
struct these results in the figure caption. However, the 
qualitative features are easily understood without any 
calculation, by folding the Fermi surface according to the 
period four charge stripe order (thin red lines in Fig. 2), 
and reconnecting the Fermi surfaces (thick black lines in 
Fig. 2) to produce gaps where ever two folded Fermi 
surfaces cross. 

As a result, so long as 2ir > q' > 37r/2, a small electron 
pocket will be formed where the Fermi surface in the 
presence of nonzero 4>n makes its closest approach to 
the (tt, 0) point. The size of this pocket is a decreasing 
function of Vq , andalsoof (q'— 3ir/2). An elclctron pocket 
with an area around 2% of the BZ (corresponding to the 
size observed in the quantum oscillations experiments) 
can be obtained by choosing small values of Vq and V2; 
for the purposes of computing quantities to compare with 
experiment, we will take take Vb = 0.02i and V2 = O.Oli 
and the remaining parameters as in Figs (a) and 2(b). 
(However, for graphical clarity, in Fig. 2, we have chosen 
larger values of Vq and V2 to make the character of the 
Fermi surface reconstruction easier to visualize.) If q is 
considerably less than 7r/2, then even for non-zero Vq (so 
long as it is not too large), there will be a hole pocket 
near the (0,7r) point in the BZ, as shown in Fig. 2(a). 
Conversely, if q is close to ir/2, the portion of the Fermi 
surface near (0, tt) is nearly nested, so even a small CDW 
potential (Vo ^ 0) will open a gap in that region of the 
BZ - the case shown in Fig. 2(b) (39|. The remaining 
portions of Fermi surface are open, as can be seen in the 
figures, and so cannot contribute to quantum oscillations 
(although they can contribute to transport anisotropy). 

From the temperature dependence of the quantum os- 
cillations, it is possible to extract the value of the effective 
mass. The effective mass that enters the expression for 
the quantum oscillations is 



2vr 



dk 



(2) 



where the integral is over the Fermi-surface of the rele- 
vant pocket, and eg is the dispersion. With this defini- 
tion, the electron pockets in Fig. 2(a) and 2(b) have effec- 
tive mass m* /m e = 0.19eV/t and 0.22eV/t, respectively, 
where m e is the bare electron mass. The measured effec- 
tive mass in the quantum oscillations is typically about 
twice the free electron mass, although it varies somewhat 



with doping concentration. This would agree with the 
above result if we set t « O.leV [io| . 

The same effective mass enters the contribution of the 
electron pocket to the specific heat. As we are ignor- 
ing all effects of bilayer coupling, there are two identi- 
cal electron pockets per Cu-0 bilayer in the reduced BZ 
(— 7r < ky < 7r and — 7r/4 < k x < 7r/4). At low temper- 
atures, this gives rise to an electron pocket contribution 
to the specific heat, C e /_ poc = 7 e ;_ poc T with "f e i- P oc = 
2.92(m*/m e ) mJ-K -mol -1 , which is comparable to the 
total specific heat measured[6[ in the relevant field range 
- for instance^, at 40T, ^ exv « 4.2 mJ-K~ 2 -mol -1 . How- 
ever, the total contribution to the specific heat from all 
the portions of the Fermi surface in the reconstructed 
band-structure is j t h = 3.28(eV/i)mJ-K _2 -mol _1 for 
Fig. 2(a) and 4.03(ey/<)mJ-KT 2 -mor 1 for Fig. 2(b), 
which for any reasonable value of t are considerably larger 
than the measured specific heat. This is a serious prob- 
lem with interpreting experiments in YBCO in terms of 
the simple model considered here. Possibly, as conjec- 
tured by the authors of Ref. 0], superconducting fluc- 
tuations still produce a residual gap on at least some 
portions of the "normal state" Fermi surface. 

The amplitude of the charge modulation, defined as 
the difference between the highest and lowest site charge 
density, is calculated to be about 0.05e in both cases of 
Fig. 2(a) and Fig. 2(b), which is reasonably close to the 
value 0.03±0.01e, adduced[I3 from the high field NMR 
experiments. The relative modulation of density of states 
at the Fermi level (p) is computed to be about 5p/p~ 0.1 
for Fig. 2(a) and Sp/p « 0.04 for Fig. 2(b), from which 
the corresponding variation of the Knight shift can be 
predicted. (The spatial average of p is p m 0.5/i.) 

Obviously there are many further effects that have not 
been taken into account in the present study. For exam- 
ple we have ignored the effects of fluctuations, and of in- 



cipient spin order. Moreover, there is some evidence [41 



from neutron scattering studies of phonon anomalies 
in YBCO, that there is a strong tendency toward 2kp 
charge density wave formation along the stripe direction 
(i.e. along the y direction), which we have neglected as 
well. Certainly, the rough consistency we have demon- 
strated with experiment is far from a unique property of 
our simple model. For example, at least as good agree- 
ment with the quantum oscillation data has been ob- 
tained under the assumption of dDW order 

In addition, there are many aspects of the this problem 
that arc perplexing, and may well imply that much more 
subtle effects are at play. A few of the more vexing are: 
1) It is not at all clear how to reconcile the existence of 
an electron pocket in the antinodal region, near (tt, 0), 
with clear spectroscopic evidence of a well developed 
pseudo-gap in this region. The NMR data suggests 
that charge stripe order condenses only when the field ex- 
ceeds a large critical value, of order 20T, and possibly this 
field is large enough to quench the antinodal pscudogap. 
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However, an intuitive link has been suggested between 
the existence of an electron pocket and the observed sign 



change of the Hall resistance and thermopower j43l |44 | 
below a critical temperature, T . This sign change can 
be observed at relatively low fields in YBCO with x 
near 1/8, where spectroscopic studies still show a well 
developed pseudo-gap [45 -47 1 . Moreover, similar signa- 
tures are seen in the Hall and thermopower 43M45I ] in 
Eu doped LSCO and LBCO, where the superconducting 
T c is suppressed and charge stripe order is stabilized by 
the LTT crystal structure, even in zero field and where, 
moreover, the existence of an antinodal pseudogap is well 
established [H]. 2) It is observed that at lower doping, 
x k, 8%, YBCO does not exhibit quantum oscillations. 
It has been suggested (Hf that this is due to the electron- 
pocket vanishing by a Lifshitz transition. While such a 
transition is possible in the present scenario, it does not 
seem natural. However, at this low doping, spin order is 
observed 4' 



- 51| , which at short distances can be thought 
of as spin-stripe order, but at long distances is character- 
ized as a spin-glass. Possibly the scattering of electron 
quasiparticles from the spin-glass order or slow fluctu- 
ations quenches the quantum oscillations. 3) Most im- 
portantly, the specific heat data(f| strongly suggest that 
substantial local superconducting order, and the corre- 
sponding gapping of the quasiparticle spectrum, survives 
even at the high fields where quantum oscillations are ob- 
served. How any of the physics we have addressed plays 
out in the presence of strong superconducting fluctua- 
tions is unclear to us at present. 4) Finally, it is im- 
portant to stress that the model we have treated omits 
many material specific details, such as effects of bilaycr 
splittings and the presence of chain bands, which even if 
they are not conceptually important, can greatly affect 
the outcomes of specific experiments. 

We thank M. -H. Julicn for allowing us to see his data 
prior to publicaiton, and M.-H. Julicn, C. Proust, J. 
Tranquada, E. Fradkin, and R.-H. He for insightful dis- 
cussions. DHL and HY were supported by DOE grant 
number DE-AC02-05CH11231 at UCB and SK was sup- 
ported, in part, by DE-FG02-06ER46287. 



[1] N. Doiron-Leyraud et al, Nature 447, 565 (2007). 

[2] D. LeBoeuf et al, Nature 450, 533 (2007). 

[3] S. E. Sebastian et al, Nature 454, 200 (2008). 

[4] B. J. Ramshaw et al., Nature Physics advanced online 
publication, (2010). 

[5] C. Jaudet et al, Phys. Rev. Lett. 1 00, 187005 (2008). 

[6] S. C. Riggs et aL larXiv:1008.1568l 

[7] S. Chakravarty, Science 319, 5864 (2008). 

[8] B. Vignolle et al, Nature 455, 952 (2008). 

[9] S. Chakravarty and H.-Y. Kee, Proc. Natl. Acad. Sci. 
U.S.A. 105, 8835 (2008). 
[10] S. E. Sebastian et al, Phys. Rev. B 81, 214524 (2010). 
[11] S. A. Kivelson et al, Rev. Mod. Phys. 75, 1201 (2003). 



[12] A. J. Millis and M. R. Norman, Phys. Rev. B 76, 220503 
(2007). 

[13] T. Wu, H. Mayaffre, S. Kramer, M. Horvatic, C. 
Berthier, W. N. Hardy, R. Liang, D. A. Bonn, and M-H. 
Julien, "Magnetic-field-induced stripe order in the high 
temperature superconductor YBa2Cu30j/' unpublished 
manuscript. 

[14] There is another form of CDW order that is consistent 
with the NMR data which has ordering vector (■k/2,tt), 
in which there is a charge modulation with period 2 along 
the stripe direction on top of an underlying period 4 site- 
centered stripe state. 

[15] J. M. Tranquada et al., Nature 375, 561 (1995). 

[16] P. Abbamonte et al, Nature Phys . 1, 155 (2005). 

[17] M. Huecker et al larXiv:1005.5T9Tl 

[18] J. Fink et al, Phys. Rev. B 79, 100502 (2009) and J. 

Fink et al, larXiv: 101 1.51011 
[19] B. Lake et al, Nature 415, 299 (2002). 
[20] Vershinin et al, Science 303, 1995 (2004). 
[21] C. V. Parker et al, Nature 468, 677 (2010). 
[22] Y. Kohsaka et al, Nature 454, 1072 (2008). 
[23] C. Howald et al, Phys. Rev. B 67, 014533 (2003). 
[24] J. E. Hoffman et al, Science 295, 466 (2002). 
[25] S. A. Kivelson, E. Fradkin, and V. J. Emery, 1998, Nature 

393, 550. 

[26] Y. Ando et al, Phys. Rev. Lett. 88, 137005 (2002). 

[27] V. Hinkov et al, Science 319, 597 (2008). 

[28] R. Daou et al, Nature 463, 519 (2010). 

[29] M. J. Lawler et al, Nature 466, 347 (2010). 

[30] M. Vojta, Adv. Phys. 58, 564 (2009). 

[31] M. Hashimoto et al, Nature Physics, 6, 414 (2010). 

[32] R.-H. He et al . New J. Phys. 1 3, 013031 (2011). 

[33] H. Yagi et aL larXiv:1002. 06551 

[34] T. Dahm et al, Nature Physics 5, 217 (2009). 

[35] I. S. Elfimov, G. A. Sawatzky, and A. Damascelli, Phys. 

Rev. B 77, 060504(R) (2008). 
[36] Y. Sassa et al. larXiv:1101. 56441 

[37] E. Fradkin et al, Annu. Rev. Condens. Matter. Phys. 1, 
153 (2010). 

[38] Y-J. Kao and H-Y. Kee, Phys. Rev. B 72, 024502 (2005). 

[39] A third, qualitatively distinct situation arises when q is 
considerably larger than n/2, but this does not appear to 
be pertinent to the cuprates, so it will not be discussed 
here. 

[40] A value t — 0.38eV has been obtained from first principle 
calculations by O. K. Andersen et al, J. Phys. Chem. 
Solids 56, 1573 (1995). But interactions are expected to 
renormalize t down to a smaller value. 

[41] D. Reznik et al, Phys. Rev. B 78, 094507 (2008). 

[42] In the supplemental material of R.-H. He et al, Science 
331, 1579 (2011), it is shown that some features of the 
pseudo-gap spectrum in Bi2201 can be caricatured by a 
simple model of the sort considered here, although in that 
paper, the effects of charge stripe or nematic order were 
considered by themselves, but not in combination. 

[43] J. Chang et al, Ph ys. Rev. Lett. 10 4, 057005 (2010). 

[44] F. Laliberte et al.. larXiv:1102. 09841 

[45] J. M. Tranquada et al, Phys. Rev. B 81, 060506 (2010). 

[46] T. Valla et al, Science 314, 1914 (2006). 

[47] R.-H. He et al, Nature Physics 5, 119 (2009). 

[48] D. LeBoeuf et al, Phys. Rev. B 83, 054056 (2011). 

[49] J. E. Sonier et al, Phys. Rev. B, 76, 064522 (2007). 

[50] S. Sanna et al, Phys. Rev. Lett. 93, 104507 (2004). 

[51] D. Haug et al, New J. Phys. 12 105006 (2010). 



